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ABSTRACT 


A variety of optimization algorithms for engineering 
synthesis are included in a new general-purpose optimization 
computer program called ADS-1l (Automated Design Synthesis, 
Version 1). Preliminary testing of all presently available 
algorithms is conducted utilizing several carefully selected 
problems of significant size and complexity. These include 
a problem with 56 design variables and over 3500 inequality 
constraints. 

The capabilities and utility of the ADS program coupled 
with a structural analysis' code utilizing finite element tech- 
niques is demonstratec ancl numerical results are presented 
that compare the relative efficiency and reliability of the 
various optimization algorithms. The number of function and 
gradient calculations are considered important measures of 
Merit in comparing the various algorithms. 

A comparison of results with another existing optimiza- 
tion computer code is included to document the accuracy and 
metiaoility of the ADS program. Preliminary testing of the 
ADS program demonstrates the flexibility a design engineer 
would have in selecting an optimization algorithm best 


Suited to solve a particular problem. 
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I. INTRODUCTION 


A. BACKGROUND 

The concept of structural synthesis, a new general ap- 
proach to structural optimization, wasS popularized by Schmit 
in 1960 [Ref. 1]. Structural synthesis, simply stated, 
couples finite element structural analysis with non-linear 
Mathematical programming techniques. Schmit reasoned that the 
design of structures for minimum weight was, after all, 
Simply the classic problem of allocation of scarce resources. 
He emphasized the importance of considering a multiplicity of 
distinct loading conditions and the need for inequality con- 
straints to deal with a variety of different failure modes 
Simultaneously, as well as side constraints (or bounds) on 
the size of the elements in the structure [Ref. 2]. 

Numerical techniques to solve the general non-linear, 
inequality constrained optimization problem developed rapidly 
after 1960. It was the advances of the high speed digital 
computer however, that allowed the science to fully mature. 
In fact the state of the art in mathematical programming is 
such, that the design engineer today should not find it nec- 
essary to develop his own computer program considering the 
widely available existing codes and the prohibitive costs of 
develcping a new optimization code. The state of the art 
in finite element analysis has also enjoyed a considerable 


advancement. Thus there exists today the ability to efficiently 


eZ 





design complex structures with many design variables under 
multiple loading conditions subject to a variety of con- 
straints including stress, displacement, buckling and fre- 
quency as examples. 

Structural synthesis continues to be the subject of 
active research; two specific areas for further study have 
been identified by Vanderplaats in [Ref. 3]. First is the 
need for pubiic availability of a computer code incorporating 
a variety of optimization algorithms that reflect the state 
of the art in optimization. Secondly, the efficiency, relia- 
bility and accuracy of the various algorithms need to be com- 
pared and the results well documented. With this information, 
the engineer who may not have written his own optimization 
code, would be able to intelligently select the appropriate 
algorithm with only a basic knowledge of structural synthe- 
sis concepts, and tailor the algorithm to suit a particular 
problem. 

The ADS library of design optimization algorithms was 
developed by Vanderplaats in response to the first need for 
a new general-purpose optimization computer code [Ref. 4]. 
ADS is unique insofar as it incorporates in a Single program, 
a variety of different optimization algorithms. The purpose 
of this research is to perform some of the preliminary testing 
of this code and document the comparative studies. The 
specific objectives of this thesis as well as the details of 
the development of the various computer codes will be dis- 


cussed in the remaining sections of this chapter. 


An 





Bb. tHESTS OBJECTIVES 

The primary objective of this thesis is to conduct the 
preliminary testing of all presently available algorithms in 
the ADS program. Although more than two years in develop- 
ment, numerous programming bugs remain to be ferreted out. 
Furthermore, some algorithms had never been tested with 
problems of significant size. Various default values for 
control parameters will also be determined by the preliminary 
testing. 

While testing is in progress a second primary objective 
is to compare and document the efficiency of programming, 
reliability of results and accuracy of solutions of the 
various algorithms. To insure validity of the comparative 
study all testing is to be accomplished in accordance with 
the following requirements: 

1. The same person is to test all algorithms on the 
same computer. The mainframe computer used in this 
research is an IBM 3033 system 370. 

2. Default values will be used in the comparative studies. 
"Fine tuning" of algorithms by overriding default 
settings will be avoided insofar as possible. 

3. Test cases of significant size and complexity will 
be selected for their potential to demonstrate the 
utility and flexibility of the ADS program and not 
because of their known ability to work well ona 


given algorithm. 
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Finally, a secondary objective is to compare results to 
the solutions provided by CONMIN [Ref. 5], a Fortran program 
for constrained minimization developed by Vanderplaats in 
1973. CONMIN is considered well tested and reliable; the 


comparison thus rendered should lend credence to the results. 


C. DEVELOPMENT OF COMPUTER PROGRAMS 
fioeovo> lL aucemarted Design Symchesis, Version 1) 

The primary motivation behind ADS-1 was the need to 
provide a selection of optimization algorithms in a sophis- 
ticated computer code that could be applied to a variety of 
design problems. The ability to easily override default 
values of control parameters further enhances the flexi- 
bility of the program to be tailored to suit the particular 
design problem at hand. 

The ADS program [Ref. 4] is written in subroutine 
form, well documented internally, and contains pseudo- 
dynamic dimensioning to maximize the efficient use of storage 
in the computer. Due to its inherent modularity the program 
1s easy to interrupt and restart and amenable to multi-level 
optimization. These features add to its portability and 
reflect the state of the art in modern programming practices. 

COPES, the control program for invoking CONMIN [Ref. 
6], was modified for use with ADS and is named "COPESA", 
whereby data transfer into and out of ADS is readily accomplished. 

The solution of an optimization problem is divided 


into three user defined levels: 


AS 





ie STRATEGY=-The method of optimization used may be 
direct, where control is transferred directly to the 
Surin zch, Ousineirect as In various penalty function 
methods. A complete list of strategies is in Table i. 

2. OPTIMIZER--Algorithms presently include methods for 
unconstrained functions as well as direct methods for 
constrained methods. A complete list of optimizers is 
in Laole sir 

3. ONE-DIMENSIONAL SEARCH--The user iS given a choice 
of curve fitting a polynomial with or without finding 
bounds, using the Golden Section method or uSing a 
combination of polynomial and Golden Section methods. 
A complete list of one-dimensional search techniques 
mo la faple Til. 

The program assumes the user is knowledgeable enough 
to select an appropriate combination of strategy, optimizer 
and one-dimensional search. For example, it would not be 
appropriate to use a variable metric cptimizer on a con- 
strained optimization problem unless one of the penalty 
function strategies was specified. Table IV lists the 
available options and feasible combinations are indicated. 

Zeal (Structural Analysis and Design--Trusses) 

The primary purpose of SADT by Fitzgerald in [Ref. 
36] was the development of a finite element code for three- 
dimensional indeterminate truss analysis and design. The 


code waS written such that it could be easily coupled to an 
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ISTRAT 


6 * 


7* 


8x 


PAB ee 


Strategy Options in ADS 


STRATEGY TO BE USED 
None. Go directly to the optimizer. 


Sequential unconstrained minimization using 
the quadratic exterior penalty function 
method [Refs. 7 and 8]. 


Sequential unconstrained minimization using 
the linear extended interior penalty function 
method [Refs. 9 through 11]. 


Sequential unconstrained minimization using 
the quadratic extended interior penalty function 
method [Ref. 12]. 


Sequential unconstrained minimization using 
the cubic extended interior penalty function 
method [Refs. 13 and 14]. 


Augmented Lagrange multiplier method 
[Refs. 15 through 19]. 


Sequential Linear Programming [Refs. 20 and 21]. 


Method of Centers (Method of Inscribed 
Hyperspheres) [Ref. 22]. 


Powell's Variable Metric Method for Constrained 
Minimization [Refs. 17, 23 and 24}. 


* 
Not available as of February, 1983 
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Optimizer Options in ADS 


OP err ok ko sre eon. 


None. Go directly to one-dimensional search. 
This option should be used only for program 
development. 


Method of Feasible Directions (MFD) for con- 
strained minimization [{Refs. 25 and 26]. 


Fletcher-Reeves algorithm for unconstrained 
Ininamazacion [Ref. 27]. 


Robust Method of Feasible Directions for con- 
strained minimization [Ref. 28]. 


Davidon-Fletcher-Powell (DFP) variable metric 
method for unconstrained minimization [Refs. 29 

and 30]. 

Broydon~Fletcher-Goldfarb~-Shanno (BFGS) variable 
metric method for unconstrained minimization [{Refs. 
Slee nmough 34 |; 

Random Search for unconstrained minimization. 
Random Search for constrained minimization. 


Newton's Method for unconstrained minimizaticn. 


Quadratic Programming [Ref. 35}. 


* 
Not available as of February, 1983 
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One-Dimensional Search Options in ADS 


ONE-DIMENSIONAL SEARCH OPTION [Refs. 7 and 52] 


Find brackets on the minimum of an unconstrained 
unmet don. 


Find the minimum of an unconstrained function 
uSing the Golden Section method. 


Find the minimum of an unconstrained function 
using the Golden Section method, followed by 
cubic polynomial interpolation. 


Find the minimum of an unconstrained function 
by first finding bounds and then using 
polynomial interpolation. 


Find the minimum of an unconstrained function 
by polynomial interpolation/extrapolation without 
Piense Tia narmgeooundas On the Solution. 


Find brackets on the minimum of a constraired 
Buln e tal On. 


Find the minimum of a constrained function 
using the Golden Section method. 


Find the minimum of a constrained function 
using the Golden Section method, followed by 
cubic polynomial interpolation. 


Find the minimum of a constrained function by 
first finding bounds and then using poiynomial 
interpolation. 


Find the minimum of a constrained function by 


polynomial interpolation/extrapolation without 
first finding bounds on the solution. 
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STRATEGY 


6 * 
7* 


Rx 


ONE-D SEARCH 


1 


2 


10 


Xx 


x 


Allowed Combination 
Combination Not Allowed 


TABLE IV 


Program Options in ADS 


OP iT ir ar in 
HE 2 3 > 
X X X X 
0 X 0 X 
0 X 0 X 
0 X 0 X 
0 XK 0 X 
0 X 0 X 
X 0 X 0 
X 0 X 0 
Xx 0 x 0 
0 0 0 0 
0 X 0 X 
0 X 0 X 
0 X 0 X 
0 X 0 X 
0 0 X 0 
X 0 X 0 
Xx 0 X 0 
X 0 X 0 
X 0 Xx 0 


6 * 


Not Available as of February, 
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optimizer for comparative studies. A secondary objective 
was to provide a user-friendly computer code that could be 
employed for truss analysis only. SADT was therefore 
selected as the analysis code for test cases involving 
trusses and space towers. 

Design variables may include member element cross 
sectional areas, nodal Res suns esee or both. A well written 
user's manual is included in [Ref. 36] and provides neces- 
sary details for coupling the program to an optimizer as 
well as for test case data preparation. 

The finite element method of analysis is used for 
Static analysis, and eigenvalues are computed according to 
the subspace iteration technique when frequency constraints 
are specified [Ref. 37]. Multiple static loading conditions 
can be accommodated as well as constraints on stress, Euler 
buckling, displacement and the first fundamental frequency 
of the structure. The objective function iS minimum weight 
of the structure. Side constraints may be imposed on the 
upper and/or lower bounds of the design variables. Design 
Wereiable linking 2S permitted for both member areas and 
coordinates. The user may specify different materials for 
the various members. All loads are assumed concentrated at 
the joints and the truss is treated as a discrete, pin- 


connected structure. 
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Bee PREVIOUS COMPARATIVE STUDIES 

Even though many methods are available for solving the 
constrained, non-linear optimization problem there has been 
relatively little research done in the way of comparative 
Studies since the inception of structural synthesis in 1960. 

Colville, in a landmark study in 1968, sent eight con- 
Strained croblems (three to 16 estan variables each) to the 
developers of 30 different codes. Solution times as well as 
preparation time and the number of function and constraint 
evaluations were requested from each participant [Ref. 38]. 
Colville placed great emphasis on solution times and there- 
fore developed a standard timing routine in an attempt to 
normalize solution times to eliminate differences among 
computers. He could not of course, eliminate the differences 
in the developers' abilities to efficiently code their 
problems for solution. 

Eason and Fenton tested 13 different problems on 20 dif- 
ferent codes in 1972 [Ref. 39]. They effectively eliminated 
the problems evident in Colville's study. All of their test 
case problems however, had fewer than seven independent design 
variables. 

Sangren and Ragsdell conducted a comparative study on 30 
problems in [Ref. 40]. The number of design variables in 
this study range from two to 48 while the number of con- 
straints range from zero to 19. 

The problems selected for comparative study in this 


research have from 5 to 56 design variables and from ll to 


oe 





3550 constraints, the largest problem being the design of a 
234-bar space tower subject to constraints on stress, Euler 


buckling, and displacement of joints. 
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A. OPTIMIZATION CONCEPTS 
The general, non~linear, constrained optimization 


problem can be stated mathematically as: 


Pome. Ze 
F(X) (28) 
Subject to: 
G(X) < 0 j = 1,NCON 2c ) 
ke = 1 ND (22) 
1- i-- il 
F(X) MomGalledeune Obyecerye funCEION. Lt zs the function 


with respect to which the design is optimized. It may be a 
impear OF non-linear function of the design variables x. 
Generally speaking, the objective function may be implicit 
Serexplicit functions of x. It 1S important however, that 
these functions be continuous and have continuous first 
derivatives in x. The 5 ON inequalities define the con- 
straints which the user imposes on the design. Equation 2.3 
defines side constraints or bounds on the design and are the 
limits over which F(X) and G(X) are defined. If the inequality 
Senale1on of equation 2.2 is not met for any constraint, 


me@at COnStraint 1S said to be violated. If the equality 
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Senadteremmer Equation 2.2 1s met then the constraint is 
called active. 

The ability to deal with equality constraints is alsc 
included in the ADS program. This feature was not fully 
operational at the time of this writing however, and there- 
fore was not tested. 

The n-dimensional space eoanned by the design variabies 
x 1s referred to as the design space. Any design satisfying 
equations 2.2 and 2.3 is a feasible design and the minimum 
feasible design is said to be optimal. Problems in optimi- 
zation may be classified according to whether or not they 
are constrained. Algorithms to solve these problems are 
therefore generally classified by the type of problem they 
were developed to solve efficiently. In the remaining sec- 
tions of this chapter the algorithms used in the preliminary 
testing of the ADS library will be discussed. Techniques to 
solve the unconstrained minimization problem will be discussed 
first, followed by constrained minimization methods. Lastiy, 
the various technigues for minimizing functions of cne varia- 
ble, the so-called one-dimensional search, will be discussed. 
These technigues are called upon by both major categories 
of algorithms to solve a sub-problem in the optimization task, 
wherein the following recursive relationship is commonly 


employed: 


xt = xT l 4 yest (2.4) 


ZS 





in this equation q is the iteration number, a* is the scalar 


step size and S is the vector search direction. 


B. UNCONSTRAINED MINIMIZATION 
een roauce 10N 

In the general case of unconstrained minimization of 
a multi-variable function, the calculus requires for a 
minimum solution, that the gradient of the objective function 
with respect to the design variables equat€ to zero and that 
the Hessian matrix of second partial derivatives of the 
objective function with respect to the design variables be 
positive definite (all eigenvalues > 0). If the Hessian 
matrix 1S poSitive definite a relative minimum at least is 
guaranteed. Unconstrained methods are therefore, intrinsically 
concerned with gradient information; as a result, they are 
classified according to the type of derivative information 
they require. JZero-order methods such as Random Search and 
Powell's Conjugate Directions Method are non-gradient methods 
whereas first-order methods such as Fletcher-Reeves require 
first derivative information only and so on. These methods 
as well as the variable metric methods of Davidon-rletcher- 
Powell and Broydon-Fletcher-Goldfarb-Shanno will be discussed 
in the next few sections. 

2. Non-Gradient Methods 
a. Random Search 
Random Search methods represent the simplest 


possible approach to optimization, wherein a randomly 
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selected large number of possible X vectors are evaluated 
for values of the objective functions. The X vector corres- 
ponding to the least objective function is the optimal design. 
There are many drawbacks, not the least of which is efficiency. 
The necessity to evaluate a large number of possible designs 
is required to insure a precise optimum has been obtained. 
The need to improve efficiency is the motivation behinc. many 
of the modifications available for random search methocs. 
These methods lend themselves well to coding on a hand-held 
calculator, furthermore they require little storage on the 
computer, making them efficient from that point of view. 

b. Powell's Conjugate Directions Methoed 

Powell's method is certainly the most popu-.ar, 
if not the most efficient, of all zero-order methods. Powell's 
Method is based on the concept of conjugate directions. The 
algorithm requires an initial search in n-orthogonal direc- 
tions wherein each search updates the X vector according to 
equation 2.4. 

The new search direction is found by simply con- 
necting the first and last design points; this becomes the 
n+l conjugate search direction. Powell's Method breaks down 
if a search direction makes no improvement because subsequent 
search directions will not be conjugate. A second well 
recognized problem is the tendency after a few iterations for 
the search directions to become nearly parallel. Powell 


offers a sophisticated technique to overcome this second 
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problem [Ref. 41]. Simply restarting the process with uni- 
directional searches is an effective, if not elegarit, way 

of dealing with this problem as noted in [Ref. 42]. Powell's 
Method is not presently available in ADS. 

The next logical step in sophistication is to 
provide gradient information to the optimizer. In the 
following sections the Fletcher-Reeves algorithm and variable 
metric methods will be discussed insofar as they are first- 
order methods presently available in ADS. 

3. Gradient Methods 
a. Fletcher-Reeves Method of Conjugate Dizxections 

The Fletcher-Reeves algorithm is actually a 
modification of the steepest descent algorithm with a signi- 
ficant improvement in the rate of convergence. Th: basic 


approach is to pick conjugate search directions according to: 


1Q 
ll 


CW 


TE pst (2.5) 


where: 


D 
ll 


IVP (XY) |*7| vr (x) 77] * (2.6) 


The initial search direction is in the direction 


of steepest descent: 


st = - vF(x?) (2am) 
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The method is conceptually similar to Poweli's 
Method, except now each search direction 1S conjugate. 
Theoretically, convergence for a quadratic function inn or 
fewer iterations can be expected, however, restarting the 
process every few iterations as in Powell's Method is 
usually required. 

b. Variable Metric Methods 

Variable Metric Methods retain information about 

previous iterations also. In these methods a !natrix Ht is 


created which approximates the inverse of the Hessian matrix. 


The search direction is defined at iteration q as follows: 


A Hoe) (oe) 


~s 
~ 


Again the initial search direction is determined 
by the method of steepest descent. At the end of iteration 
eee ne H Matrix 1s updated according to: 


~s 


qatt 2 yi pd (249) 


~s 


where D? is a symmetric matrix determined according to the 


~— 


PELOowing formulation: 


ih 2 2.10 
D [o + 9t/o lop? ( ) 


ww me 


the terms in this equation are defined as: 
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~~ ~s 


y = vwr(x4) - vr(x27*) ei) 
o = pry (2.13) 
1 = y Hly (2.14) 


and §@ iS a parameter used to select the form of the update 
formula, equation 2.10. The Davidon-Fletcher-Powell Method 
sets @ = 0 in equation 2.10 whereas the Broydon-Fletcher- 
Goldfarb-Shanno Method sets 6 = 1 [Ref. 42]. There are other 
possible algorithms in the class of variable metric methods 
but these two methods are the most popula and are presently 


available in ADS. 


C. CONSTRAINED MINIMIZATION 
i, Wighepatole bogota Kola! 

Constrained methods of minimization were developed 
to deal with problems that have limitations placed on a set 
of functions of the design variables. These limitations may 
be side constraints which directly impose bounds on the de- 
Sign variables, or so-called behavior constraints which are 
functions of the design variables. Behavior constraints may 
take the form of equality or inequality constraints, but in 
either case the design must satisfy the behavior constraints 


while staying within the bounds imposed by the side constraints. 
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Direct methods consider the constraints as limiting hyper- 
surfaces and attempt to directly minimize the objective 
function in their presence. In contrast, the so-called 
penalty function methods transform the constrained minimiza- 
tion problem into a sequence of unconstrained minimization 
problems. Although direct methods are often more efficient, 
indirect methods are popular because they are simple to 
invoke. The engineer must employ an appropriate unconstrained 
minimization algorithm when using a penalty function method. 

The indirect methods utilizing penalty function 
techniques may be further classified into two broad categories: 
interior and exterior. Interior methods are designed to ap- 
proach the optimum from the feasible region whereas the ex- 
terior methods approach the solution from the infeasible 
sector. A pseudo-objective function is created by imposing 
a penalty for violated constraints. The general technique 
is to minimize this pseudo-objective function as an uncon- 
strained problem. The methods require repetitive solution 
to a series of unconstrained problems thus the term, "Sequen- 
tial Unconstrained Minimization Techniques" (SUMT), is applied 
to this broad class of indirect methods. 

2. Direct Methods 

Most optimization algorithms proceed iteratively 
toward a solution from a user supplied initial xX vector 
which may or may not define a feasible design. The design 


is modified according to the recursive relationship: 
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X = x14 gts (Zs) 


where g 1s the iteration number, S 1S a vector search direc- 
tion in the design space and the scalar, a*, defines the 
distance the optimizer moves in the search direction S. 
The choice of S 1s such that the objective function is re- 
duced. The efficiency and reliability of a given optimization 
algorithm is largely due to the fundamental method of deter- 
mination of the search direction S and the step size a*. 
These methods will be discussed in the next few sections of 
this chapter. 
a. Method of Feasible Directions 

Optimization in the Method of Feasible Directions 
proceeds in two basic steps, first a usable-feasible search 
direction is determined, then a one-dimensional search 1s 
performed in this direction to reduce the objective as much 
as possible without violating constraints. The method as- 
Sumes that the initial X vector of design variables defines 
a feasible design. A usable-feasible search direction to 
improve this design is found by solving the following sub- 


problem: 


Maximize: 8 C226)) 
Subject to: 


VF(X)"S +B < 0 (2.17) 


a2 





VG(X)-S + 0,-8 < 0 Jed (2.18) 


eo) i (2.19) 


~~ 


where J is the set of currently active constraints, a = 0. 
V is the gradient operator and the components of 6 are re- 
mmrea EG 2S puSsh-off factors, which act to push the design 
away from currently active constraints. A value of unity 

for §@ will yield a search direction which approximately bi- 
sects the usable-feasible sector. 

If the initial design is infeasible it is possible 
to find a search direction that will direct the design to the 
feasible region [Ref. 42}. 

Using equation 2.19 with equations 2.16 through 
2.18 results ina linear problem of finding S except for one 
G@uadratic constraint. Zoutendijk in [Ref. 25] provides a 
mieect approach to overcome t:his difficulty. A detailed 
explanation of these techniques is provided in [Ref. 42]. 

The methcd then proceeds to update the design in 
accordance with equation 2.15. This step 1s commonly per- 
formed by polynomial interpolation but a variety of one- 
dimensional search methods may be used. 

b. Robust Method of Feasible Directions 

The Robust Method of Feasible Directions is a 
new algorithm presently being developed by Vanderplaats, and 
incorporates the best features of the Method of Feasible 


Directions (MFD) and the Generalized Reduced Gradient (GRG) 
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Method [Ref. 28]. Only gradients of active constraints are 
required in the MFD, which is considered an attractive feature, 
while the GRG method has the nice feature of precisely follow- 
ing the constraint boundaries from one vertex to the next 
without the need to move away from the constraints. The 

Robust MFD retains these desirable features but does not re- 
quire the addition of slack variables peculiar to the GRG 
method, thus avoiding the large matrix operations associated 
with the GRG method. The method involves solving the following 


search direction sub-problem: 


Maximize: 


- VE(X)°S (Ze 
SUbJECE CO: 

VG.(X)*S <0 ed (2221) 

2S (2522) 


This is the same form as the direction finding 
sub-problem in MFD except the dimensionality is reduced by 
the elimination of the variable 8. The following advantages 
in determining the search direction in this manner are re- 
peated here from [Ref. 28] for convenience: 

1. The dimensionality of the design problem is not in- 
creased by the addition of slack variables to the 


inequality constraints. 
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Ze wee algorithm fom finding S 1s specifically designed 
for inequality constrained problems, thus improving 
ein) CLeEnCYy;: 

3. Only gradients of active constraints are required. 

4. The number of dependent variables is greatly reduced 
in comparison to the GRG method, thus a reduction in 
the size of the sub-problem in the one-dimensional 
search is achieved. 

Equality constraints are effectively handled as 
a special case of inequality constraints. Initially infeasi- 
ble designs require a modification to the search direction- 
finding sub-problem where the violated constraints are treated 
as inequality constraints. A direction to the feasible region 
is then determined in a manner similar to the Method of 
Feasible Directions. 

The Robust method incorporates a particularly 
attractive feature of infrequent gradient calculations. 
That is, gradients of active constraints are treated as con- 
stants for several iterations thus greatly reducing the 
computational cost of the algorithm. It should be noted that 
1f infrequent gradient calculations are not used the method 
yields the same results as the GRG Method. 

The one-dimensional search is performed in the 
same manner as for the GRG method. Significant in this 
procedure is the fact that Newton's Method is employed to 


drive the active constraints corresponding to the dependent 
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variables to zero. This procedure usually requires several 
iterations. 
The Robust Method of Feasible Directions shares 
some of the limitations of the GRG method [Ref. 28]. 
l. It produces infeasible designs and relies on Newton's 
Method to return to the feasible region. 
woe has difficulty dealing with highly non-linear 
BUMeELONS . 
3. If the analysis is itself iterative the method may 
be unable to satisfy constraints due to the resulting 
mastabii 1 tye 
Sa  undirect Methods 
ADS incorporates several SUMT methods, namely, 
exterior, extended interior, and Augmented Lagrange Multi- 
plier (ALM) penalty function methods. The numerical ill- 
conditioning often encountered in SUMT methods is reduced in 
the ALM method. This method has therefore received wide 
attention in the literature and is included in the ADS library. 
Ali SUMT methods create a pseudo-objective functicn 


of the general form: 


OX, = F(X) + r_P(X (2223) 
(2 p) (X) p x) 
where F(X) is the original objective function, P(X) is the 


penalty function and the multiplier, eo? determines the 


Magnitude of the penalty applied. The following sections 
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discuss in more detail the technique of determining P(X) 
which is the fundamental basis of each method. 
a. Exterior Penalty Function Method 
The basic mathematical formula for determining 


the penalty function P(X) is: 


a 2 elle 2 
P(X) = ) {MAX[0,q.(X)]}“ + J) thy (X) ] (2724) 
: = a k=l * ~ 


A penalty is imposed if, and only if, an inequality, G(X), 
Sarequality, Hy (X), Cemstralnt 12S yiolated. The "offending" 
constraint is squared to provide a slope of zero for the 
penalty function at the constraint boundary thus insuring a 
continuous first derivative for the pseudo-objective func- 
tion. The second derivative is not required to be contin- 
uous hewever, therefore if second-order methods are employed 
PMmene UnNCcCOnstrained minimization, numerical ill-conditioning 
Inmay result [Ref. 42]. 

The multiplier, ee is critical in this method as it 
is in all SUMT methods [If _ 1s chosen small the pseudo- 
objective function is easily minimized but may result in 
extreme constraint violation; whereas a large a will guard 
against this, the resulting problem is usually numerically 
ill-conditioned. Therefore the algorithm starts with a small 
ro which is then increased by a factor y. At each iteration 


¢6 is minimized starting from the previous optimum solution. 
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As r 1s increased in the sequential optimiza- 
tion process, the pseudo~-objective function becomes increas- 
ingly non-linear. The constrained optimum solution is also 
approached from the infeasible region. In other words the 
optimum is approached with a series of infeasible designs, 
none of which are usable. The interior penalty function 
method approaches the optimum from the feasible sector with 
a series of improving feasible designs. This attractive 
feature is discussed in the next section. 

b. Interior Penalty Function Method 
The most common formulation for the penalty 


meneccign in this method is: 


P(X) = 


Iie39g 


[-1./g, (x) ] (2.25) 


j=1 


resulting in a more complicated pseudo-objective function to 


minimize: 


- 
Reh) = F(X) + EPO] te Y the] (2.26) 


'O 
nw 


Note that equality constraints (hy) are dealt 
with in the same manner by interior and exterior methods. 
The significant difference between the methods, besides the 
menmmulation of P(X), is the fact that in interior methods 


the penalty parameter, a5 is sequentially decreased with 
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every SUMT iteration, while in exterior methods Es 1s sequen- 
tially increased. Interior methods result in the approach of 
the optimum solution from the feasible region as es 2 Op but 
is discontinuous at constraint boundaries. The exterior 
method, on the other hand, is well-defined everywhere, but 
leads to an optimum solution only in the limit as a ror, 
The extended interior penalty methods are designed to incor- 
porate the best features of both methods by effecting a 
transition between the interior and exterior methods at a 
point in the optimization task. Needless to say, this transi- 
tion point 1s critical and therefore of fundamental concern 
in the various extended interior penalty function methods, 
which are discussed next. 
c. Extended Interior Penalty Function Method 

The chief advantage of the interior penalty method 
is that it results in a sequence of improving feasible de- 
ens £rOm an initially acceptable starting point. This 
desirable feature is maintained in this method by a judicious 
selection of the parameter, ¢«, in the formulation of the 


Bemalty function P(X): 


mM 
PC) VS es 227) 
~ j=1 3) ne 
where 
ora = rae ee lg CS ne (2.28) 
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9, (X) = “2. - 9, (x) I/e° alg Bib) > We (277 29)) 


The parameter, ¢€, is a small negative number and 
Signifies the transition from the interior to the exterior 
methods [Ref. 42]. These equations define the linear ex- 
tended interior penalty function. Because the second deriva- 
tive of ao is discontinuous, Haftka and Starnes 


created the quadratic extended interior penalty function by 


Smanging equation 2.29 to: 


g(x) = -1./e{{g,(x)/el® - 3.19, (X)/e] + 3.3 (2.30) 


las Get = %e 
J ~ 
Again the degree of non-linearity of ¢ is in- 
creased as a price for the second-order continuity. 
The linear and quadratic extended interior 
penalty methods are both critically dependent on the selec- 
tion of ¢«. Haftka and Starnes recommend that « be determined 


Secording to: 
6 = “C (08) © 1/3 ae) 2 (2.31) 


where C is a constant. At the beginning e« is chosen in the 
mamge —.3 < € < -.1 and is chosen such that the objective 
and pseudo-objective functions are equal; the resultant 


value of C is thus determined [Ref. 43]. 
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The quadratic extended interior penalty method 
has the disadvantage that the penalty increases dramatically 
for badly violated constraints. The variable penalty function 
method attempts to overcome this difficulty while continuing 
to insure second order continuity at the transition point. 
The selection of € in the variable penalty method is recom- 


mended by Prasad in [Ref. 44} as follows: 
q 
Se = = roe ee 
B ( 5 ( ) 
where 
Wazrowme <— Gg < 1/S Lor Sa eo (2753) 


and 8 iS a poSitive constant chosen such that ¢« is initially 
near zero. In ADS, the variable penalty method is used 
wherein S = 3 thus the strategy is referred to as the cubic 
extended interior penalty function method. 
d. Augmented Lagrange Multiplier Method 

The efficiency of SUMT methods can be improved 
by the inclusion of Lagrange multipliers, thus reducing 
dependency of the algorithm on the choice of the penalty 
parameters. The Lagrangian is created for equality con- 


strained problems as follows: 


Z 
oe) eee) 
Sa , k=1 


A, A (X) (2.34) 
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Since the minimum of the Lagrangian provides the 
solution to the general equality constrained problem, a 
pseudo-objective function, called the augmented Lagrangian 


is created using the exterior penalty function method: 
- 2 
A(X,A,r) = F(X) + eK + xh, (YT) Cee) 


The method starts with the following values for 


= +1. if vh,(X)-VF(X) < 0 (2.36) 
= -l. if Vh,(X)-VF(X) > 0 (2.37) 


The pseudo-objective function, A(X,AyLp), ls then 
minimized holding ro and ,} constant. A new set of Lagrange 


multipliers is calculated according to: 


ee 2x hy (xP) ok = Lye (2.38) 
The parameter i is sequentially increased as in 
the exterior SUMT method and the unconstrained minimization 
problem is solved for ry and }. The process is repeated 
until convergence is achieved. 
The method is easily extended to handle inequality 


constraints by converting them to equivalent equality 





constraints by the addition of slack variables. A more com- 
plicated augmented Lagrangian is then formulated as the 


pseudo-objective function: 


m 
2 2 2 
Pak ,A,2,r_) = F(X) + ie (Gexy eae) ero a (Xk) FS 
eS” p 3 re i ee oe 
(2?) 
where there are m slack variables, ae These are calculated 


as a sub-problem and so do not increase the dimensionality 
of the optimization task. Note that the pseudo-objective 
function has continuous first derivatives with respect to x 
but discontinuous second derivatives at a WS = aa thus 
second order techniques should be avoided in the unconstrained 
minimization problem. The method has several attractive 
features repeated here from [Ref. 42]. 
1. The method is relatively insensitive to oe accordingly 
it 1s not necessary to increase XD COee: 
2. Equality constraints and inequality constraints pre- 
cisely equal to zero are possible. 
3. Acceleration to an optimum is achieved by updating 
the Lagrange multipliers. 
4. The starting point may be feasible or infeasible. 
5. At the optimum any Lagrange multiplier not equal 


eemzer@ Will identify an active constraint. 
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4, Other Constrained Minimization Methods 

The discussion of optimization algorithms has been 
restricted to non-linear programming techniques insofar as 
these methods are fully operational in ADS. Sequential 
Linear Programming is another category of optimization tech- 
niques which will be included in the ADS library where a 
particular problem is linearized and a solution sought for 
the resulting linear approximation. Considering these tech- 
niques are, in theory, well-developed and quite effective 
eos additional capability will enhance the utility of ADS. 

The basic approach is to linearize the objective and 
Gemetraint functions and obtain a solution to this approxima- 
tion using the algorithm developed for linear programming. 
The process is iterative and therefore the techniques are 
referred to as Sequential Linear Programming (SLP). It is 
pointed out in [{Ref. 42] that fully constrained problems 
usually converge rapidly while under-constrained problems 
often have difficulty in converging to an optimum solution. 
The difficulty may be overcome somewhat by sequential reduc- 
tion of move limits on the optimizer. SLP characteristically 
produces a sequence of improving infeasible designs. The 
Method of Centers also called Method of Inscribed Hyperspheres, 
has the dual advantage of approaching the optimum with a se- 
quence of improving feasible designs while following a path 
down the "center" of the design space. This method is dis- 


emesocad in the remainder of this section. 
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The basic approach in the Method of Centers is to 
inscribe a hypersphere in n-dimensional design space created 
when all of the constraints and objective function are 
linearized. The design then moves to the center of the 
hypersphere. This procedure is repeated to convergence 
within some user-specified tolerance. In the case of under- 
constrained problems the method is subject to the same problem 


as SLP in imposing move limits on the optimizer. 


BD. FUNCTIONS OF ONE VARIABLE: THE ONE-DIMENSIONAL SEARCH 
Po citcoduct lon 

The one-dimensional search, as it is commonly referred 
Pemin algorithms for optimization, usually applies to deter- 
mining a*, the step size to be taken in the search direction 
S. Finding the minimum of any function of one variable is 
Simply finding the point at which the first derivative 
vanishes. Since the function is not always an easily obtained 
analytic function in optimization, it is necessary to make 
some fundamental assumptions so that appropriate numerical 
analysis techniques may be brought to bear. Accordingly, the 
functions are assumed unimodal, that is, the function has 
only one relative minima in the region of concern. The func- 
tions are also assumed continuous as are their first and 
second derivatives. These assumptions will assure convergence 
fond minimum. 

In the remaining sections of this chapter the methods 


used to conduct the one-dimensional search will be discussed. 
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2. Polynomial Approximation 


The basic procedure in the polynomial approximation 
method is to evaluate the function at several points and 
then fit a polynomial curve to the data points using an 
appropriate curve-fitting technique. The minimum of this 
curve is approximately equal to the minimum of the true 
Function. The method is simple, requires only a few func~ 
tion evaluations and is generally reliable for functions 
which are not too highly non-linear. 

It is well known that a higher order polynomial will 
Fit the data points more accurately; this gain in accuracy 
however can complicate the process of finding the minimum of 
the resulting polynomial. Also, interpolation between points 
1s preferred to extrapolation beyond the region enclosed by 
the data points. The process of finding the minimum of the 
polynomial requires finding the point where the first deriva- 
tive vanishes. Alternatively, there are numerical analysis 
techniques available to find the minimum or zeros of a higher 
order polynomial. These methods are not discussed here. 

3. Golden Section Method 

The Golden Section Method is popular because the rate 
of convergence is Known and the requirements for function 
unimodality and continuity are relaxes. The disadvantage 
of the method lies in the inherently large number of function 
evaluations required as compared to other one-dimensional 


search methods. 
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The method involves picking two intermediate points, 


x) and X51 between given upper and lower bounds, Xo and Xe 
such that xy < X, - The function is then evaluated at Ay and 


xX. and one of the previous bounds is replaced by one of the 


intermediate points as follows: 


F(X, ) > F(X.) Xo eS (2 40) 
+ ( 
F(X.) > F(X)) x Xo meh Ll) 
The process is repeated until some user specified 
tolerance is satisfied. Fundamental to the method is the 
selection of the intermediate points. The Golden Section 
number, 1.61803, is used for this purpose: 
(X. - Xo) /(Xy - X,) =, 1261303 fee42) 


The Golden Section provides the ideal sequence for 
dividing the interval such that the minimum number of function 
evaluations is required. The advantage of this method is 
guaranteed accuracy whereas the relatively large number of 
function values required is a distinct disadvantage. 

A similar method, the Fibonacci Search, based on the 
series of Fibonacci numbers, traps the minimum in successively 
smaller intervals. The Fibonacci Search is occasionally 


more efficient than Golden Section but is far more complicated. 
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4. Finding Bounds on the Solution 
This method is usually used to obtain brackets on 
the solution, then Golden Section or polynomiai methods are 
called to complete the one-dimensional search. 
The method begins with an assumed initial lower 


bound X, and a proposed upper bound xX These two points are 


ie 


then evaluated as F(X, ) and F(X). lige F(X) > BCX, then i 
is the true upper bound. Assuming the slope of the function 


at Xp is negative at X the solution is complete. If F(X) 


(a 
= F(X, ) then the following update formula is applied itera- 
tively to achieve the desired bounds: 

eld Old 


new _ _ : 
xy = skew ict a) xX) a Xo (23-415) 


new old (2.44) 


where a = golden section number = 1.61803. 

Note that if the last three values of this iterative 
procedure are retained along with the function values, the 
three required points by the Golden Section and Polynomial 
methods are already available. 

Many algorithms (e.g., MFD) require the constrained 
minimum of F(X). Polynomial and Golden Section methods are 
also used in ADS for this purpose. Note that the X used here 


1s actually a* in equation 2.4. 


48 





Mine nee INARY TESTING OF ADS—1 


eee 6=6LNTRODUCTION 

Selection of test problems in a comparative study is of 
primary importance. Considering one of the objectives of 
this thesis is to demonstrate the utility and flexibility of 
the ADS library, test cases were selected from two fundamen- 
tally different areas in which optimization is commcnly 
used. These areas are structural design (trusses, frames, 
Space towers, etc.), and ship synthesis. There are many 
other areas in engineering where optimization is employed 
but the areas chosen here are selected for comparative study 
in this research due to the availability of the analysis codes. 

A good test case 1S one in which no single const:raint 
dominates the design. Three different truss cases were 
selected that met this criteria. They were also chosen 
because they are significant in size and complexity and thus 
would demonstrate the comparative efficiency and reliability 
of the various algorithms to be tested. Truss cases are 
popular in the literature because differences due to modeling 
details and idealizations can be eliminated easily; they also 
lend themselves well to finite element methods of analysis. 
Because analytically optimum solutions to the test cases 
are not available, solutions obtained by the well developed 
and thoroughly tested optimization program "CONMIN" are 


provided as a base-line for the results from ADS. 





The remaining two test problems consist of a 10-variable 
cantilever beam optimized for minimum volume, and the con- 
ceptual design of the FFG-7 Perry Class Frigate where the 
objective function is taken to be the full load displacement 
of the vessel. The detailsof the varicus test cases are 


presented in the following sections. 


gee DESCRIPTION OF TEST PROBLEMS 
1. 10-Variable Cantilevered Beam 

The 1l0O-variable cantilever beam test case was devel- 
oped by Vanderplaats in 1979 as a teaching aid for a graduate 
level course in DeSign Optimization. The problem is quite 
simple, yet the solution is not easily obtained. The beam 
consists of a specified number of equal length sections; 
each section has a rectangular cross section with the height 
constrained not to exceed 20 times the width. This equates 
to a crude buckling constraint. The maximum stress at the 


left end of each section is constrained as follows: 
Oo. oe Osa i Wes (2 eelr) 


The beam is cantilevered and tip loaded with a force 
of 10 kips downward and the total tip deflection is con- 
strained not to exceed two inches. Material properties of 
the beam conform to steel where Young's modulus, E = 30 x 20° 
psl. 

The initial X vector of design variables consisting 


of height and width dimensions of each section is tabulated 
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in Appendix A. There are five equal length segments in the 
overall length of 200 inches, resulting in an initial volume 
of 8000 in. The objective function is the minimum volume 
subject to the constraints of stress, displacement, and 
height to width ratio. A three dimensional drawing of the 
initial design is shown in Figure 3.1 and the optimum solution 
is shown in Figure 3.2. The tip deflection of the optimum 
beam is actually two inches downward but no attempt is made 
to show this in the figure. 
2. 10-Bar Planar Truss 

Numerous test cases for planar trusses (2-dimensional) 
and space towers (3-dimensional) can be found in the litera- 
ture. In particular, the 10-bar planar truss has been used 
in [Ref. 45] to demonstrate how the stress-ratio method, 
which seeks a fully-stressed design, yields poor results when 
members with significantly different allowable stresses are 
specified [Ref. 46]. 

The configuration of the 10-bar planar cantilever 
Eruss is shown in Figure 3.3 and is subject to a single load 
condition of 100 kips downward at nodes two and four. The 
initial cross-sectional areas of the truss elements and 
bounds on the areas are listed in Appendix B. 

There are 20 constraints consisting of maximum and 


minimum stresses in each of the 10 members as follows: 


eect Coy as 625 ksi i = 1-8,10 (3.2) 
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Ste kcl agen < to0) ksi oe (323) 


where 1 is the memper element number. It should be noted 
that element nine in Figure 3.3 has twice the allowable 
stress of the other members. The objective function is 
Minimum weight of the structure. Material properties include 
Young's Modulus, E = naeeenO. psi and y = .l lb/in? corres- 
Bonding to the properties of aluminum. 

3. Conceptual Design FFG-7 PERRY Class Frigate 

The details of this test case may be obtained in 
[Ref. 47] where Jenkins optimized the conceptual design of 
a FFG-7 Perry Class Frigate. More specifically, he coupled 
Ene Reed synthesis model for surface combatant ships [Ref. 
48], with the non-linear optimize CONMIN, a FORTRAN program 
for constrained function minimization, via the control pro- 
gram COPES. COPES/CONMIN was developed in 1973 by Vanderplaats 
[Ref. 5] and has been used in a variety of engineering appli- 
cations. The objective function is the full load displace- 
ment of the vessel. 

The design variables used in the preliminary testing 
of the ADS program are the same as those used by Jenkins: 
accordingly a comparison of results is appropriate. The 
independent design variables include: 

1. LBP - Length between perpendiculars, ft. 
2. L/B - Length to beam ratio 


3. B/H - Beam to draft ratio 
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4. Cp - Prismatic coefficient 
5. Cx - Midship section coefficient 

The initial values of these variables as well as 
their upper and lower bounds are listed in Appendix C. 

There are 13 constraints on the design, these are explained 
in detail in [Ref. 47], and are not repeated here. 
4. 47-Bar Planar Tower 

The 47-Bar planar tower shown in Figure 3.4 was 
introduced in the literature in [{Ref. 49], wherein the tower 
was designed subject to multiple loading conditions. The 
Same tower was designed for optimum geometry in [Ref. 50] 
subject to stress and Euler buckling. In [Ref. 50] sub- 
structuring was also used. The two sub-structures were over- 
lapped so that several members were in both sub-structures. 
Meer. 21) presents configuration optimization with the addition 
of frequency constraints. 

The 47-Bar planar tower used in this research 1s 
discussed in the remainder of this section. Initial cross- 
sectional areas of the truss elements, nodal coordinates and 
bounds on these parameters are tabulated in Appendix D as 
well as the details regarding displacement constraints and 
loading conditions. Steel was selected as the material for 
all members with Young's Modulus, E = 30 «10° psi and 
i= .3 lb/in>. 

All elements are subject to the following constraints 


on stress: 
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where i is the element number. Tubular members are speci- 
mheed with a Glameter Eo thickness ({(d/t) ratio = 10. Euler 
buckling is prohibited by constraining the buckling stress 


in the members according to the following equation: 


co. > oO = -10.1nEA,/8L‘ t= ea es) 
Finally, the first fundamental frequency of the 
Structure is required to exceed 5. cps. Two non-structural 


weights of 500 lbs each are attached at nodes 17 and 22 to 


rh 


acilitate the eigenvalue problem solution. 

Member areas and coordinates are linked to maintain 
symmetry about the vertical Y axis. Nodes 15, 16, 17 and 
22 are fixed in space and nodes 1 and 2 are constrained to 
meron the X axis. The resulting problem thus reduces to 27 
member sizing variables and 17 configuration variables for 
a total of 44 independent design variables and 436 constraints 
on stress, Euler buckling, displacement and frequency. 

The optimum design is shown in Figure 3.5. It should 
be noted there was no attempt to show member sizing variables 
mieene figures. 

5. 234-Bar Space Tower 
The configuration of the 234-Bar space tower is shown 


in Figure 3.6. Initial cross-sectional areas of the truss 
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Figure 3.5. Optimum Design at the 47-Bar Planar Tower 
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elements, noaal coordinates, bounds on these parameters, as 
well as the details regarding displacement constraints and 
loading conditions are tabulated in Appendix E. Aluminum 
was selected as the material for all members with Young's 
Modulus, E = 10 x 10° psamanciey = .1 lb/in?. 

All elements are subject to the following constraints 


on stress: 
Soensil  < Os eco ksa 1 = €,234 (32/0) 


where 1 is the element number. Tubular members are speci- 
fied with a diameter to thickness (d/t) ratio = 10. Euler 
buckling is prohibited by constraining the buckling stress 


in the members according to the following equation: 
é — 
C > Cy = ol LO.17EA, /8L. ii — ae (Sa) 


Member areas and coordinates are linked to maintain 
Symmetry about the vertical Y axis. Nodes l, 2, 3 and 4 are 
constrained to lie on the XZ plane. The resulting problem 
thus reduces to 56 member sizing variables and 3550 constraints 


on stress, Euler buckling, and displacement. 


C. COUPLING ANALYSIS AND OPTIMIZATION COMPUTER CODES 
The test case data files were prepared in accordance 
with the user's manual for SADT [Ref. 36] and the user's 


manual for COPESA, similar to [Ref. 6]. The problems were 
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then coupled to the ADS library of optimization algorithms 
Via a brief driver program in the case of trusses and towers 
and via COPESA on the cantilever beam and Ship design cases. 
All test case results were printed and filed for future refer- 
ence. The default values for all program control parameters 
such as convergence tolerances were used insofar as possible. 
Gradients were calculated analytically for the truss and 
tower cases and by finite differences in the cases of the 
beam and ship. 

The results obtained were carefully tabulated and opti- 
mum solutions determined based on the best objective function 
and the fewest equivalent function evaluations. This param- 


eter was computed as follows: 


NFE = IFCALL + NDV*IGCALL (Sec) 


where iFCALIL is the number of objective and constraint 
function evaluations, IGCALL is the number of times gradients 
are evaluated by the user and NDV is the number of design 
Variables. This provides an equivalent number of function 
evaluations that would be required if all gradients were 
Calculated by finite differences. If gradients are calcu- 
lated by finite differences, IGCALL will be zero because 
IFCALL includes the function evaluations needed to calcu- 


late gradients. 
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IV. RESULTS AND CONCLUSIONS 


A. INTRODUCTION 

There are presently 85 possible, meaningful combinations 
ef strategy, optimizer and one-dimensional search methods 
available in the ADS library. Testing all methods on all 
problems is not practical considering some test cases con- 
sume over 40 minutes of CPU time per run. Accordingly, 
the scope of research was limited to testing all strategies 
and all optimizers with three one-dimensional searches on 
all five problems, for a total of 260 test case computer runs. 
The two one-dimensional search methods not tested were bounds 
only and polynomial without bounds. The results are tabulated 
in Appendix F. In Tables V through IX the best optimum de- 
Signs to each of the five problems are presented. Optimum 
design A represents the best objective function achieved, 
whereas optimum design B represents the solution within 5% 
of the objective function for optimum design A but which had 
the fewest equivalent function evaluations. Both solutions 


were required to have no violated constraints (g(X) < 0.01). 


B. RELATIVE RANKING OF OPTIMIZATION METHODS 
fe panccution Tame 
A timing routine available in the Non-IMSL library 
at the computer center was utilized to record execution time 


in CPU seconds for each test run. Times were then averaged 
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for all runs using the same one-dimensional search on a given 
problem. In other words, CPU time per function evaluation was 
averaged for all runs recorded on any given table in Appendix 
F. These run times, when multiplied by the equivalent number 
of function evaluations, iS a good approximation of CPU seconds 
to optimize a problem with any given combination of strategy, 
optimizer and one-dimensional search. For example, average 
CPU time per function evaluation for the cantilever beam range 
from .002581 seconds to .0037283, whereas the range on the 234- 
Bar space tower is .32508 to .36011 seconds. It is readily 
apparent that on problems of significant size, like the 234- 
Bar tower (Table XXXV, Appendix F) run times of 34 CPU minutes 
may be realized. The significant point is that the efficiency 
Siean atgorithm to reduce NFE to a minimum is of vital concern 
on problems of practical interest. 
2. Number of Function/Gradient Calculations 

A perusal of all results in Appendix F reveals that 
Girect methods are far more efficient than indirect methods 
as far as NFE is concerned when solving constrained minimi- 
zation problems. Furthermore it is apparent that the ALM 
method is effective in reducing NFE for SUMT methods as theory 
would suggest. 

Contrary to expectations, there is no apparent trend 
that would indicate which unconstrained minimization method 
is "best" to use when employing a SUMT method for the solu- 
tion of a constrained problem. Perhaps more extensive testing 


would result in establishing these desirable guidelines. 
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A review of Table VII points out an interesting fact 
concerning NFE. In this table the optimum solutions for the 
Piemmeceosencase are recorded. Note that the Method of Feasi- 
ble Directions results in a quite acceptable objective func- 
tion in 55 function evaluations while a SUMT method (exterior 
penalty) required 555 function evaluations to achieve a 
Slightly better result! This situation is shown graphically 
in Figure 4.1. The point here is for the user to be aware 
of the possibility that an optimizer may be using an inordinate 
amount of computer resources to achieve an insignificant gain 
in the objective. 

ee vatues Of the Obsyective Function 

PeGeGnpablson of @Gbyective Functions points out that 
in general all presently available algorithms are working 
well in ADS with the exception of SUMT methods on the 234- 

Bar space tower. In this case the optimizers were unable to 
overcome the constraint violations and make progress toward 

a solution; the trouble is attributed to needed refinement in 
choosing the penalty parameters. 

The efficiency, reliability, and accuracy of the 
various algorithms however, is clearly demonstrated on the 
other four test problems as recorded in Tables XXIII through 
XXXIV in Appendix F. In these four test cases, extremely 
good objective functions were obtained and generally resulted 
in the production of feasible designs (no violated constraints). 

Tables V through IX record the best objective func- 


tion achieved for each problem. Again, the direct methods 
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TABLE V 


Optimum Design of 10-Variable Cantilever Beam 


DESIGN 
VARIABLES 
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TABLE VI 


Optimum Design of 10-Bar Planar Truss 


OBJ: 


INITIAL 
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TABLE VII 


Optimum Conceptual Design of FFG-7 Perry Class Frigate 


DESIGN 
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Optimum Design of 47-Bar Planar Tower 


OPTIMUM 
DESIGN B 
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no violated 


constraints 





DESIGN 
VARIABLES 


Si 
2 
53 
34 
2 
36 
Se 
38 
0 
40 
4l 
42 
43 
44 


OBd: 


LAB ey LiL 


NTA 
VALUES 


.60000E+02 
.24000E+03 
-60000E+02 
.36000E+03 
.30000E+02 
.42000E+03 
.30000E+02 
.48000E+03 
.30000E+02 
.54000E+03 
.90000E+02 
.60000E+03 
. 30000E+02 
.60000E+03 


~-28650E+04 


te eA: 
IGCALL: 
NEE: 
TS LRAw 
POE: 
TIONED: 


ga 


(Contd) 


OPTIMUM 
DESIGN A 


~-40992E+02 
~2/7822E+0 3 
oe oO 
.37854E+03 
~24938E+0 2 
(4516288203 
~24265E+02 
-47202E+0 3 
, 2 s09 OF U2 
~51176E+03 
.82690E+02 
~5919 3E+03 
(2650 7b4-0 2 
-DoI7ZLE+O3 


- O04 aE O4 


948 
nl 
SoZ 
2 

4 

2 


OPTIMUM 
DESIGN B 





DESIGN 
VARIABLES 


Oo Oxy DW MN Sf& WH KN FF 


Ce NON NNT NN nh Ae ee ee Ee ee Be ese 
rF oN WON HD UN SF WwW NY PO YO ON HD TO F&F WD NY FP OC 


UJ 
NO 


TABLE IX 


INITIAL 
VALUES 
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yielded a better result than the indirect methods on the two 
smaller problems while the SUMT methods prevaiied on the 
larger problems. This suggests that indirect methods deal 
with a multitude of active constraints more effectively than 


direct methods. 


Ge COMPARISON WITH CONMIN 

All test cases were run on CONMIN to provide a base-line 
for the comparative studies. Results from CONMIN are recorded 
on the tables of test results in Appendix F. It is inter- 
esting, if not surprising, that with the excepticn of the 
234-bar tower test case, ADS routines were able to achieve 
better solutions than CONMIN. 

CONMIN basically utilizes a Feasible Directions ézlgorithm 
for constrained problems. The fact that a different combina- 
tion surpassed CONMIN on each test case supports the notion 
that the optimization algorithm employed should suit the prob- 
lem at hand to gain maximum efficiency. In the past the 
thrust has been to merely alter the program parameters of 
the same algorithm to deal with fundamentally different prob- 
lems. ADS now offers a convenient method for selecting an 
algorithm best suited to the problem at hand. This flexi- 
bility further enhances an engineer's ability to apply opti- 


mization concepts to the various disciplines in design. 


D. ADDITIONAL CONCLUSIONS 
Preliminary testing of the ADS library resulted in the 


modifications of several default values for the various 


74 





optimizers that improved their efficiency dramatically. As 
ADS is fully implemented additional testing will be required 
to insure all algorithms are as efficient, reliakle and 
accurate as possible. 

A difficulty with a program of this broad capability is 
to provide the user with a concise set of guidelines identi- 
fying which method or class of methods shculd be selected 
for a given problem. The problem is exacerbated by the 
selection of default values, in other words, a default value 
which may work well on one problem may cause premature con- 
vergence on a different problem. Accordingly, judicious 
selection of default values in ADS requires cons:.derable 
effort supported by extensive testing on a variei:y of problems 
as the algorithms become operational. 

Results given in Appendix F and the optimum solutions 
tabulated in Tables V through IX are an indication of relia- 
bility, to be sure, however the results are preliminary and 
the algorithms are constantly being revised and improved. 
The equivalent number of function evaluations (NFE) provide 
a measure of relative efficiency of the optimizer to achieve 
an optimum solution, the goal being to minimize the use of 
computer resources while maximizing the reduction of the 
objective function. It should be noted and is evident in 
the results tabulated that the efficiency and reliability 
are problem dependent. Therefore a wise selection of the 


appropriate algorithm and tailoring the program parameters 


15 





to suit the problem at hand is required. ADS achieves this 
flexibility and enhances the design engineer's ability to 


use optimization as a viable design tcol. 
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APPENDIX A 


1O-VARIABLE CANTILEVER BEAM TEST CASE 


Table X describes the initial X vector of independent 
design variables; the first five variables are segment 
widths and the remaining five are segment heights; side 


constraints or bounds on the variables are also included. 


TABLE X 


Initial X Vector for the 10-Variable Cantilever Beam 


SEGMENT DIMENSIONS (INCHES) 


DESIGN LOWER i he UPPER 
VARIABLE BOUND VALUE BOUND 
il - S0000E+00 e2ooo0 +O 1 SOU Or--O 
2 POUCOOETCO PAu 000R+GL PS O.O UO to OL 
S P20000E +00 we OO Ol PSOOO00ETOL 
7 . S0000E+00 Poouoor+oL PoN00E+0L 
2 .50000E+00 ZU OUR +0 1 oO GOO. OL 
6 SOK Ore .20000E+02 . LOOOOE+03 
q . LOOOOE+02 OOO aa s9) 2 . LGOO0OE+O3 
8 PEO COG EO PawOOOR+FOZ PLOQOORA Os 
9 eo OOR+O 2 PaO O00E FOZ . LODQOE+0 3 
10 PLOUOOn. OZ TOCOOE+OZ PLOOCORtO. 
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APPENDIX B 


1lO-BAR PLANAR TRUSS TEST CASE 


Table XI describes the initial X vector of independent 
design variables; the ten variables consist of the truss 
element cross-sectional areas. Side constraints or bounds 
on the variables are also included in the table. Table XII 


lists the nodal coordinates in inches. 


TABLE XI 
Initial X Vector for the 10-Bar Planar Truss 
DESIGN TRUSS 


VARIABLE ELEMENT 
NUMBER NUMBER 


GRGss SEeGCIPTONAL AREAS (SO. IN;) 


LOWER BCUNDS INITIAL VALUES UPPER BOUNDS 


i. i OO OC E00 OOO E02 - LOOQOE+04 
2 2 » LOOQDOE+00 J600CE 02 . LOOOOE+04 
= 3 s LO URIG Eo LeooeE 02 - LOOOOE+04 
3 4 » LOOOOQE+00 POO Ore 2 . LOOOOE+04 
> 5 POS UE +O 0 CUD CE+O2 . LOOOOE+04 
6 6 - LOQOQQE+00 SOO OE+ 02 - LOQOQOE+04 
A y .L0000E+00 OOOO ET OZ - LOOOOE+04 
8 8 peo 0 OOo 70.0 008-62 . LOOOOE+04 
9 3 PlOCOOE LOU PEO e Ciao 2 - LOOQOE+04 
10 10 ee OO s00 PUOe OBO 2 -lOOCCE ECs 
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Initial Nodal Coordinates of the 10-Bar 


NODE 
NUMBER 


X 


. /Z2000E+03 
a 200 On +0 3 
. 36000E+0 3 


7 oO 00R+0 3 


0 


0 


TABLE 


XITL 


Planar Truss 


COORDINATES (INCHES) 


BG 


. 36000E+03 


0 


- JG000E+O3 


0 


. 36000E+03 
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APPENDIX C 


CONCEPTUAL DESIGN OF THE FFG~-/ PERRY CLASS FRIGATE 


Table XIII lists the initial X vector of independent 
design variables; side constraints or bounds on the variables 


are included in the table. 


TAMEE ALI 


Initial X Vector for the FFG-7 Preliminary Design 


DESIGN 

VARIABLE PARAMETER LOWER INITIAL UPPER 

NUMBER BOUNDS VALUES BOUNDS 
a LBP . 30000E+03 . 30000E+03 . /O000E+03 
2 L/B . 70000E+00 .90700E+00 .L2000E+00 
3 B/H .20000E+01 .31400E+01 -40000E+01 
4 Cp .50000E+00 .59000E+00 .I0000E+00 
5 Cx< .75000E+00 . 75000E+00 .90000E+00 
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APPENDIX D 


47-BAR PLANAR TOWER TEST CASE 


Tables XIV and XV describe the initial X vector of inde- 
pendent design variables. The 27 variables in Table XIV are 
the initial elemeni: cross-sectional areas. Table XV lists 
tne initial nodal coordinates; 1/7 of these are independent 
design variables. It should be noted that since symmetry 
about the Y-axis exists only nodes on the positive side are 
listed. Side constraints or bounds on the variables are 
included in the tables. Tables XVI and XVII describe the 


loading conditions and displacement constraints respectively. 
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TABLE XIV 


Initial Member Areas for the 47-Bar Planar Tower 


DESIGN 
VARIABLE 
NUMBER 


Oo ON HD ONO & WY NY FF 


NN NN ee se ee ee 
“HW MW £& WwW NY FO HWM WON KH UH fh Ww NY F CO 
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NM PW R BR WH NN 
ow on fF NY OO © 
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LOO OOE=0\5 
eeOOCUE—O0> 
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MEO OOE—05 
LOCC r—0> 
LOC OE—0> 
Meow OO = Os 
ROUOOCE—0'S 
prOOTCE—O> 
-LO000GE-05 
PAOCuur Os 
 LOOGGE-05 
Oe Or — OS 
muOvogE=05 
- LOOOOE-05 
PLOOgOE—O> 
- LO0008-05 
EO OMOL— OS 
FLOodoE—O > 
POCO OE—0> 
PeOwCvGE =O > 
PROG UOE=O0> 
EO 0 OE —0lS 
PEO OO OE —O> 
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INITIAL 
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- LOOOQOOE+04 
. LOOOOE+04 
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TABLE XV 


Initial Nodal Coordinates of the 47-Bar Planar Tower 


DESIGN 


VARIABLE NUMBER 


NUMBER 
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COORDINATES 
(INCHES) 

xX vs 
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0 eZ OOO 3 
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TAB iS eA vt 


Loading Conditions on the 47-Bar Planar Tower 


LOAD NODE LOADS APPLIED (LBS) 
CONDITION NUMBER 14 lane 
a iby .6000E+04 -.14000E+05 
Ze 0 0 
Z i 0 0 
ae .6000E+04 =. 14000E+05 
3 ey . 6000E+04 =o our+0>5 
ee .6000E+04 -.14000E+05 
TABLE XVII 


ro 


Displacement Constraints on the 47-Bar Planar Tower 


moOAD NODE DISPLACEMENT CONSTRAINTS (INCHES) 


SOVpiTION NUMBER Diy sel ron LOWER BOUNDS 


1 iy - 5000E+01 
oe Ne - .5000E+01 
2 ie x -.5000E+01 
7 Ms - .5000E+01 
5 7 xX - . 5000E+01 
Lz Y - 5000E+01 
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APPENDIX E 


234-BAR SPACE TOWER TEST CASE 


Table XVIII describes the initial X vector of indepen- 
dent design variables which consist of the member cross 
sectional areas. Side constraints or bounds on the variables 
are included in Table XVIII. Table XIX lists the nodal 
coorcinates. Tables XX and XXI describe the loading conditions 


and cisplacement constraints respectively. 
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TABLE XVIII 


Initial Member Areas for the 234-Bar Space Tower 


DESIGN TRUSS 
VARIABLE ELEMENT 
NUMBER NUMBER 
iL 1-4 
2 ao 
3 Done 
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DESIGN TRUSS 
VARIABLE ELEMENT 
NUMBER NUMBER 
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Nodal Coordinates of the 234-Bar Space Tower 
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.84000E+02 
77 S000 FOZ 
27 ,000E+F02 
H DO C0E+O2 
Poe OO Bae 
.66000E+02 
.66000E+02 
.66000E+02 
.66000E+02 
.57000E+02 
.57000E+02 





NODE 
NUMBER 


eel 
SZ 
33 
34 
35 
36 
oy 
38 
a9 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
56 
SL 
2 
53 
54 
59 
56 
>/ 
58 


TABijp i. CCOnt 'd) 


Xx 


wo 7 OO0E702 
5 5) NOOO 0) 22 
~-48000E+02 
~-48000E+02 
.48000E+02 
-48000E+02 
Seer 02 
me OO UE OZ 
Pee WOU 0 2 
= SOMONE) 2 
PoOOOUR+ Oc 
SoU 000+ OZ 
. 30000E+02 
neUUOOE-+ OZ 
. 30000E+02 
POOUDOE FOZ 
SS UCCOr 02 
TUCCOETO2 
-90000E+02 
PICOCCETOZ 
yevOocor+02 
. 60000E+02 
-60000E+02 
soOUOU0E--0 2 
-60000E+02 
-60000E+02 
.60000E+02 
-60000E+02 


COORDINATES (INCHES) 


x 


- 84000E+03 
. 84000E+03 
.96000E+0 3 
-96000E+03 
.96000E+03 
. 96000E+0 3 
. LO800E+04 
. LO800E+04 
. LO800E+04 
. LO800E+04 
. L2Z000E+04 
~ L2Z000E+04 
. L2Z000E+04 
~L2000E+04 
.12480E+04 
.12480E+04 
~12480E+04 
~12480E+04 
~12480E+04 
.12480E+04 
.12240E+04 
.12240E+04 
~12480E+04 
.12480E+04 
~12240E+04 
.12240E+04 
.12480E+04 
. 124 80E+04 


89 


Z 


7 CUCE+ OZ 
po VOOR UZ 
-48000E+02 
-48000E+02 
-48000E+02 
-48000E+02 
sss OUOR- OZ 
rol O0OERF UZ 
22 CUCE-TOZ 
2. JUCOR 02 
, 000GR LOZ 
PO UCOORa Oe 
. 30000E+02 
. 30000E+02 
sO C000 EO2 
SOUR U2 
SJOUDCE +1 0i2 
Pe OO Cm Oz 


0) 
0 


Pio O0r 402 
SL o000B702 
LS UUCE LOZ 
oC OE OZ 
wEs100 Ob 02 
20 00E4 0 2 
J 000ETUZ 
oO 0002 





TABLE XX 


Loading Conditions on the 234-Bar Space Tower 


LOAD NODE LOADS APPLIED (LBS) 
CONDITION NUMBER FX FY FZ 
1 49 .6000E+04 -.20000E+05 
50 .6000E+04 - ,.20000E+05 
2 49 .6000E+04 - .20000E+05 
50 -.6000E+04 -.20000E+05 
3 49 .6000E+04 -.20000E+05 0 
50 . 3000E+04 -.10000E+05 ~50000E+04 
4 49 . 3000E+04 -. 10000E+05 -.50000E+04 
50 . 3000E+04 - .10000E+05 .50000E+04 
5 49 -.3000E+04 . LOOOOE+05 .50000E+04 
50 -,3000E+04 . LOO00E+05 -.50000E+04 
TABLE  XxX1 


Displacement Constraints on the 234-Bar Space Tower 


LOAD NODE DISPLACEMENT CONSTRAINTS (INCHES) 
CONDITION NUMBER Direc Dron LOWER BOUNDS UPPER BOUNDS 
L 49 AGe =. J000E+01 S00GEFOL 

50 a, Y = S00 OHO . 5000E+01 
2 49 GG = OOO a OL - SJ000E+01 
30 a, ¥ -.5000E+01 TSO O0E EOL 
3 49 Ee =. SOUR Oi sI000EFO2 
50 aed -.5000E+01 UO OU: 
49 Earp -.5000E+01 FOOUUER Us 
30 ee =e 00 ORO. wo OGL G:l 
> a ere eo 0 Oise Ou SOO E OL 
50 a, Y Se OO Ora Ol PaO VUE sO 


90 





aeeeNDek F 


ADS-1 PRELIMINARY TEST RESULTS 


PiemscsULbeomon tie Preliminary Eesting Of the algorithms 


available in version 1 of the ADS library in February 1983 


meee summarized in Tables XXIII through XXXVII. The nomen- 


clature used in these tables is defined in Table XXII. 


TABLE AA0 1 


Definition of Terms in Test Results 


TERM 


NDV 
NCON 
OBJ 
NAC 
NFE 
SUMT 


ear 
Hove sAT INT 
QUAD EXT INT 


Sp rc EXT 
INT 


ALM 
MFD 
DFP 
BFGS 


DESCRIPTION 


Number of 
Number of 
Obj ective 
Number of 
Number of 


independent design variables 
constraints on design 
function value 

active constraints 


Equivalent Function Evaluations 


Sequential Unconstrained Minimization 
Technique 


Exterior Penalty Method 


Linear Extended Interior Penalty Method 


Ovadratic 


Extended Interior Penalty Metheod 


Cubic Extended Interior Penalty Method 


Augmented Lagrange Multipliers Method 


Method of Feasible Directions 


Davidon-Fletcher-Powell Algorithm 


Broydon-Fletcher-Goldfarb-Shanno Algorithm 


opr 





TABLE XALTIT 


Test Results 10-Variable Cantilever Beam (IONED: 2,7) 


ONE-DIMENSIONAL SEARCH: Golden Section 
NDV = 10 NCON = ll (Stress, Displacement and H/B Ratio) 
AVERAGE CPU TIME PER FUNCTION EVALUATION: .27097E~-02 seconds 


OPTIMIZER: fl 2 3 4 5 
METHOD OF FLETCHER~ ROBUST D.F.P. Bb eG.e. 
STRATEGY : FEAS. DIR. REEVES Men. 
0 OBJ .47264E+04 .38513E+04 
DIRECT NAC 1 10 
NFE 594 EeZ 
il OBJ .39742E+04 .39188E+04 .39067E+04 
SUMT NAC 1* 2* 4 
(EXT)  NFE 568 691 738 
2 
SUMT OBJ .47597E+04 .39268E+04 .39591E+04 
(LIN- NAC 3 Z 3 
EXT- NFE 736 a7 795 
INT) 
3 
SUMT OBJ .4759 7E+04 .39288E+04 .39591E+-04 
(QUAD- NAC 3 il 3 
EXT- NEE 736 ae 795 
INT) 
4 
SUMT OBJ 4759 7E+04 ~39241F+04 .39591E+04 
(CUBIC NAC cS 3 3 
EXT- NEE 736 892 795 
INT) 
5 
SUMT OBJ . 3879 3E+04 .38968E+04 .38961E+04 
(ALM) NAC 4x Oe oe 
NFE 922 692 799 


CONMIN RESULTS: 
OBJ: .40808E+04 


NAC: 4 
NFE: 313 


* = VIOLATED CONSTRAINT (S) 





Test Results 10-Variable Cantilever Beam (IONED: 


TABLE XXIV 


oie) 


ONE-DIMENSIONAL SEARCH: Golden Section + Cubic Polynomial 


NDV = 10 


NCON = 


i 


METHOD OF 


FEAS. DIR. 


-49 338E+04 
1 
403 


-40808E+04 


OPTIMIZER: 
STRATEGY : 

0 OBJ 

DIRECT NAC 
NFE 

i. OBJ 
SUMI NAC 
(EXT) NFE 

2 
SUMT OBJ 
(LIN- NAC 
EXT- MME 
INT) 

3 
SUMI OBJ 
(QUAD- NAC 
EXT- NEE 
INT) 

4 
SUMI CBI 
(CUBIC NAC 
EXT- NFE 
INT) 

5 OBJ 
SUMT NAC 
(ALM) NFE 
CONMIN RESULTS: 
OBJ: 

NAC : 4 
NFE: 313 


FLETCHER- 


. 39 856E+04 
1s 
440 


-4739 7E+04 
5 
549 


-47394E+04 
S 
aoe 


-47395E+04 
3 
a2 


. 39 740E+04 
2* 
770 


* 


S 
ROBUST 
Mer a; 


. 38513E+04 
10 
365 


-25810E-02 Seconds 


= 
Deer. 


-40222E+04 
5 
611 


- 39563E+04 
3 
594 


. 3956 7E+04 
3 
599 


. 59566704 
3 
599 


- 3865 7Et04 
Q* 
587 


ll (Stress, Displacement and H/B Ratio) 
AVERAGE CPU TIME PER FUNCTION EVALUATION: 


Bee Ges. 


. 389 84E+104 
1% 
501 


- 3920 5E+04 
i 
ogg 


eso Zio 


L 
599 


og 2 ers 


90 


- 38654E+04 


4* 
766 


VIOLATED CONSTRAINT (S) 





Test Results 10-Variable Cantilever Beam (IONED: 


TABLE XXV 


ONE-DIMENSIONAL SEARCH: Bounds + Polynomial 


NDV = 10 


OPTIMIZER 
STRATEGY : 
0 OBJ 
DIRECT NAC 
NFE 
1 OBJ 
SUMT NAC 
(EXT)  NFE 
2 
SUMT CBJ 
(LIN- NAC 
EXT- \NFE 
INT) 
5 
SUMT OBI 
(QUAD- NAC 
EXT- FE 
INT) 
4 
SUMT CBJ 
(CUBIC NAC 
EXT- NFE 
INT) 
5 OJB 
SUMT NAC 
(ALM) NFE 


NCON = 


I 


METHOD OF 
FEAS. DIR. 


-49772E+04 
i 
287 


CONMIN RESULTS: 


OBJ: -40808E+04 
NAC: 4 
NEE: 313 


ll (Stress, 
AVERAGE CPU TIME PER FUNCTION EVALUATION: 


BEE ene k= 


. 39589E+04 
2 
Sy) 7) 


-50299E+04 
i 
342 


. 49 392E+04 
uE 
400 


-49982E+04 
1 
369 


. 39678E+04 
+ 
Dol 


* =- 


. 4 
ROBUST Dawe « 
per. OD. 


. 3851 7E+04 
10 
Le 


. 39489E+04 
1* 
333 


. 39495E+04 
2 
490 


Oo orn04 
2 
388 


. 39260E+04 
2 
436 


. 38536E+04 
2* 
396 


94 


a9} 


Displacement and H/B Ratio) 
»~37283E-02 Seconds 


Geo - 


. 39849 E+0 4 
Q* 
295 


pO om 04 


2 


Sabo 


. 39689E+04 


434 


. 39562E+C4 
a: 


41] 


. 38580E+04 
3% 
461 


VIOLATED CONSTRAINT (S) 





TABLE AXVI 


Test Results 10-Bar Planar Truss (IONED: 2,7) 


ONE-DIMENSIONAL SEARCH: Golden Section 
HDvy = 10 NCON = 20 (Stress) 


AVERAGE CPU TIME PER FUNCTION EVALUATION: .41632E-02 Seconds 


OPTIMIZER: il Z 5 4 5 
METHOD OF £=FLETCHER- FOBUST DEES, BeenaGeo. 

STRATEGY: FEAS. DIR. REEVES Mote D. 

0 OBJ .15436E+04 - L4955E+04 
DIRECT NAC 1 8 

NFE 677 319 

il OBJ . 1670 8E+04 ~15558E+04 .15272E+04 
SUMT NAC 6 5 8 
(EXT) NEE 860 824 967 

2 
SUMT OBJ . 16205E+04 ~L5104E+04 .15750E+04 
(LIN- NAC 5 6 5 

EXT- NFE ES S29 15) 

INT) 

3 
SUMT OBJ -16205E+04 -L5967E+HC4 .15750E+04 
(QUAD- NAC 3 4 5 

EXI- NFE ee, L356 eo 52), 

INT) 

4 
SUMT CBI . L6205E+04 ~L5967E+04 .15750E+04 
(CUBIC NAC 3 4 5 

EXT- NFE teES7 1356 Jeoae 

INT) 

5 OBJ ~15733E+04 ~L5285E+04 .14987E+04 
SUMT NAC y 8 8 
(ALM) NEE 1138 1091 1389 
CONMIN RESULTS: <— VIOLATED sCONSTRAINT (S) 
OB : -L5009E+04 
NAC: 10 
NFE: G14 





TABLE XXVIII 


Test Results 1O-Bar Planar Truss (IONED: ?,8) 
ONE-DIMENSIONAL SEARCH: Golden Section + Cubic Polynomial 
NDV = 10 NCON = 20 (Stress) 


AVERAGE CPU TIME PER FUNCTION EVALUATION: .393298-02 Seconds 


OPTIMIZER ils 2 3 4 5 
METHOD OF FLETCHER- ROBUST Dat cP, BiakeiG.S:. 

STRATEGY: FEAS. DER. REEVES Mere). 

0 OBJ .15445E+04 . 14960E+04 
DIRECT NAC 1 8 

NFE 666 244 

1 OBJ .16717E+04 ~15743E+04 .15390E+04 
SUMT NAC 6 6 7 
(EXT) NFE 622 721 681 

Z 
SUMT OBJ . 16160E+0 4 -15102E+04 .16357E+04 
(LIN- NAC 3 5 2 

EXT- NFE 885 105% 856 

INT) 

3 
SUMT' OBJ . 16202E+04 -1576GE+04 .15684E+04 
(QUAD- NAC 3 5 5 

EXT=- NFE 826 1002 902 

INT) 

4 
SUMT OBI . 161 85E+04 ~L5747E+H04 .16352E+04 
(CUBIC NAC 3 5 4 

EXT- NFE 936 37 1016 

INT) 

5 OBJ - 160 76E+04 -14996E+04 .149965+04 
SUMT NAC 7 8 8 
(ALM) NFE 786 868 812 


CONMIN RESULTS: 


OB °: 
NAC: 
NEE: 


.15009E+04 
10 
414 


x — 


96 


VIOLATED CONSTRAINT (S) 





TABLE XXVIIT 


Test Results 10-Bar Planar Truss 


ONE-DIMENSIONAL SEARCH: 
20 (Stress) 
AVERAGE CPU TIME PER FUNCTION EVALUATION: 


NDV = 10 


OPTIMIZER 
STRATEGY: 
0 OBJ 
DIRECT NAC 
NFE 
i OBJ 
SUMT NAC 
(EXT)  NFE 
2 
SUML OBJ 
(LIN- NAC 
EXT- QMFE 
INT) 
3 
SUMT CBJ 
(QUAD- NAC 
EXT- NEE 
INT) 
4 
SUMT OBJ 
(CUBIC NAC 
EXT- MFE 
INT) 
5 OBJ 
SUMI NAC 
(AIM) NFE 


NCON = 


i 


METHOD OF 


FEAS. DIR. 


» L5313E+04 
sl 
471 


CONMIN RESULTS: 


OBJ: -L5009E+04 
NAC: 10 
NFE: 414 


PEeIGHER= 


. 16488E+04 
6 
384 


. 16034E+04 
5 
676 


. L6809E+04 
~ 
525 


. L6822E+04 
2 
BvD 


.15265E+04 
Ax 
423 


xk — 


3 
ROBUST 
een. 


. 14974E+04 
8 
2a 


oF 


(TONED: 


Bounds + Poiynomial 


-43674E-02 Seconds 


DIS ee 


L67SaE+04 
“ 
S10 


. LolO2E+04 


DN tad 


on 


.L5113E+04 
4 
o2il 


.1510-4E+04 
4 
a2 


-15414E+04 
6* 
446 


= VIOLATED CONSTRAINT (S) 


Aa, D)) 


Ben. Geol 


~15752E+04 
8 
Soo 


.15095E+04 
- 
807 


- L5LOOE+04 
S 
831 


. L5008E+04 
8 
UZ 


vee O4 
8 
780 





TABLE XXI4 


Test Results Conceptual Design FFG-7 (IONED: 2,7) 


ONE-DIMENSIONAL SEARCH: Golden Section 
NDV = 5 NGON = 13 


AVERAGE CPU TIME PER FUNCTION EVALUATION: 2257 /4E—-0)) Seconds 


OPTIMIZER: 1 2 3 4 5 
METHOD OF FLETCHER- ROBUST ODF... B.F.G.S 
STRATEGY:  FEAS. DIR. REEVES M.F.D. 
0 OBJ .35122E+04 .35077E+04 
DIRECT NAC 1 2 
NFE 155 170 

i OBJ .35114E+04 .35044E+04 .35039E+04 
sumr NAC 2 2 

(EXT) NFE 51 514 555 

2 
SUMT OBJ .35224E+04 .35060E+04 .35067E+04 
(LIN- NAC i i i 

EXT- NEE 502 544 604 

INT) 

3 
SUMT OBJ .35265E+04 .35087E+04 .35052E+04 
(QUAD- NAC 1 1 i 

EXT- FE 504 537 597 

INT) 

4 
SuMT = CBJ .35260E+04 .35122E+04 .35088E+04 
(CUBIC NAC 1 1 1 

EXT- FE 503 514 535 

INT) 

5 oBJ .35107E+04 .35063E+04 .35082E+04 
SUMT NAC i O* 2 

(ALM) NFE 285 284 248 


CONMIN RESULTS: 


OBJ: 
NAC: 
NFE: 


60 


oe Ob O4 
3 


* = VIOLATED CONSTRAINT (S) 





TABLE XXX 


Test Results Conceptual Design FFG-7 (IONED: 3,8) 
ONE-DIMENSIONAL SEARCH: Golden Section + Cubic Polynomial 
NDV = 5 NCON = 13 


AVERAGE CPU TIME PER FUNCTION EVALUATION: .299550E-01 Seconds 


OPTIMIZER: i 2 3 4 5 
METHOD OF FLETCHER- kOBUS'T Dee. B.F.G.S. 
STRATEGY: FEAS. DIR. REEVES M.F.D. 
0 OBJ .35112E+04 .35079E+04 
DIRECT NAC 0O 2 
NFE 117 107 
il OBIT .35L1LOE+04 .35039E+04 .35043E+04 
SUMI NAC a 2 2 
(EXT) NFE 374 388 393 
DE 
SUMT OBJ .35235E+04 .35088E+04 .35089F+04 
(LIN- NAC i A i 
EXT- NFE 359 366 368 
INT) 
3 
SUMT OBJ .35206E+04 .35088E+04 .35091F+04 
(QUAD- NAC i al 1 
EXT- NEE 352 356 351 
INT) 
4 
SUMI CBI .35261E+04 .35082F+04 .35087E+04 
(CUBIC NAC a 1 1 
EXT- NFE 354 399 358 
INT) 
= OBJ .35130E+04 ~35065E+04 .35073E+04 
SUMT NAC Q* 1* 2 
(ALM) NEE 192 204 210 


CONMIN RESULTS: 


OBJ : 
NAC: 
NFE: 


60 


~SO 266404 
3 


* — 


= Vio DD MGONSTRRAINT (S) 





Test Results Conceptual Design FFG-7 


TABLE XXXI 


ONE-DIMENSIONAL SEARCH: Bounds + Polynomial 


NDV = 5 


OPTIMIZER: 
STRATEGY: 

0 CBI 
DIRECT NAC 

NFE 

ih OBJ 
SUMI NAC 
(EXT)  NFE 

2 
SUMT OBJ 
(LIN- NAC 
EXT- NFE 
INT) 

3 
SUMT OBJ 
(QUAD- NAC 
EXT- NFE 
INT) 

4 
SUMT OBJ 
(CUBIC NAC 
EXT- NFE 
INT) 

5 OBJ 
SUMI NAC 
(ALM) NFE 


NCON = 13 
AVERAGE CPUT TIME PER FUNCTION EVALUATION: 


iL 
METHOD OF 
FEAS. DIR. 


. 35120E+04 
i 
55 


CONMIN RESULTS: 


GBo : 
NAC: 3 
NFE: 60 


PoolZSE+04 


REEVES 


oOo mrO4 


2 


soe bor. 04 
0 
238 


. 35305G+04 
1 
227 


~ 39.135E104 
i 
228 


. 35086E+04 
Ox 
130 


eo 


3 
ROBUST 
Mon.) 


. 350 78E+04 
2 
78 


100 


Done. 


. 350 85E+04 
Z 
262 


. 35174E+04 
al 
242 


. 35084E+04 
Hi 
220 


. 35090E+04 
1 
228 


. 35085E+04 
1 
Lee 


= VIOLATED CONSTRAINT (S) 


(IONED: 


ee) 


-28926E-01 Seconds 


tue or 


- 35083E+04 
i 
260 


» 35150E+04 
a 
242 


. 35091 E+04 
df 
226 


. 3509 0E+04 
1 
222 


. 35095E+04 
a 
174 





BABE E exxkT T 


Test Results 47-Bar Planar Tower (IONED: 2,7) 
ONE-DIMENSIONAL SEARCH: Golden Section 

NDV = 44 
AVERAGE CPU TIME PER FUNCTION EVALUATION: 


NCON = 436 (Stress, Displacement, Buckling, and Frequency. 
.53264E-0l Seconds 


OPTIMIZER: a 2 3 4 5 
METHOD OF FLETCHER: ROBUST Dee. B.F.G.S. 
STRATEGY : FEAS. DIR. REEVES M.F.D. 
0 OBJ .60646E+04 .40012E+04 
DIRECT NAC 0 10 
NFE 756 1443 

il CBI .59466E+(4 .36788E+04 .35326E+04 
SUMT NAC 2 9 ila 

(EXT) NEE 947 2154 2348 

2. 
SUMT OBJ 479 85E+04 ~30141E+04 .33997E+04 
(LIN- NAC 3 12 10 

EXT- NFE 2819 3192 2618 

INT) 

3 
SUMT OBJ .53544E+04 .30141E+04 .33997E+04 
(QUAD- NAC 4 2 10 

EXT- NFE 2207 3192 2618 

INT) 

4 
SUMT OBJ .53544E+04 .30141E+04 .33997E+04 
(CUBIC NAC 4 ie 10 

EXT~ NFE 2a 3192 2618 

INT) 

5 OBJ .40389E+04 .23645E+04 .23832E+04 
SUMI NAC o* Q* 13* 

(AIM) NFE Waa 7260 5434 


CONMIN RESULTS: * = VIOLATED CONSTRAINT (S) 


OBJ: .38078E+04 
NAC: 9 
NEE: ZO? 1. 


wok 





TABLE XXXIII 


Test Results 47-Bar Planar Tower (IONED: 3,8) 

ONE-DIMENSIONAL SEARCH: Golden Section + Cubic Polynomial 
NDV = 44 NCON = 436 
AVERAGE CPU TIME 


(Stress, Displacement, Buckling, and Frequency) 
PER FUNCTION EVALUATION: ~47939E-Ol1l Seconds 


OPTIMIZER: E 2 3 4 5 
METHOD OF FLETCHER= FOBUST D.F.P B.F2G.S. 

STRATEGY: FEAS. DIR. REEVES Mir =D. 

0 OBJ .61268E+C4 .45382E+04 
DIRECT NAC 2 Z 

NFE 404 749 

au OBJ .55562E+04 .36696E+04 .36743E+04 
SUMT NAC 3 8 5 
(EXT) NFE 1458 2390 2020 

2 
SUMT OBJ -46477E+04 -36015E+04 .35362E+04 
(LIN- NAC 4 oh 10 

EXT= NFE 24€5 2450 2408 

INT) 

3 
SUMT OBJ .44886E+04 .36063E+04 .35387E+04 
(QUAD- NAC 5 a 8 

EXT= NEE 2438 ya53 224] 

INT) 

4 
SUMT OBJ .45720E+04 .36102E+04 .35357E+04 
(CUBIC NAC 3 6 10 

EXT= NFE 2064 2252 2404 

INT) 

5 OBJ .9 /343E+04 .32050E+04 .26490E+04 
SUMT NAC 2 5 5* 
(ALM) NFE 883 3816 6258 


CONMIN RESULTS: 


OBJ: .38078E+04 
NAC: 9 
NFE: 202 1 


ee 


eZ 


= VIOLATED CONSTRAINT (S) 





TABLE XXXIV 


Test Results 47-Bar Planar Tower (IONED: 4,9) 
ONE-DIMENSIONAL SEARCH: Bounds + Polynomial 
NDV = 44 NCON = 436 


AVERAGE CPU TIM& PER FUNCTION EVALUATION: 


(Stress, Displacement, Buckling, 


~-43639E-01 Seconds 


and Frequency) 


OPTIMIZER: L 2 3 = 5 
METHOD OF FLETCHER- ROBUST BrP. BatniGeo. 

STRATEGY : FEAS. DIR. REEVES Monae. 

0 OBJ .611i15E+04 -44852E+04 
DIRECT NAC 2 2 

NEE 523 970 

a OBJ .61499E+04 .50367E+04 .59280E+04 
SUMT NAC 2 4 Z 
(EXT) NFE oS LLOL 726 

2 
SUMT OBI .62486E+04 .97233E+04 .58948E+04 
(LIN- NAC 2 0 Z 

EXT- NFE 1214 1059 1174 

INT) 

3 
SUMT OBJ .62344E+04 .36405E+04 .53255E+04 
(QUAD- NAC 0 4 Z 

Pa NEE m3 2600 1314 

INT) 

a 
SUM? OBI .09038E+04 .58969E+04 .562918+04 
(CUBIC NAC 2 Z Z 

Eat NEE 1741 EZ 1226 

INT) 

> OBI .61200E+04 .24724E+04 *** 
SUMT NAC 4 4x cs 
(ALM) NFE 871 3459 ae 


CONMIN RESULTS: 


OBJ: .38078E+04 
NAC: 9 
NFE: Ae 


*% 


VIOLATED CONSTRAINT (S) 


Ors 





TABLE XXXV 


Test Results 234-Bar Space Tower (IONED: 2,7) 


ONE-DIMENSIONAL SEARCH: Golden Section 
NDV = 56 NCON = 3550 (Stress, Displacement, and Buckling) 
AVERAGE CPU TIME PER FUNCTION EVALUATION: . 3601LLE+00 Seconds 


OPTIMIZER: ili 2 3 4 5 
METHOD OF FLETCHER- ROBUST Deer. B.F.G.S. 
STRATEGY: FEAS. DIR. REEVES M.F.D. 
) OBJ .75378E+05 .53526E+05 
DIRECT NAC 1 2 
NFE 297 1545 
il OBJ .48933E+05 .52357E+05 .46937E+05 
SUMI NAC 1* Q* O* 
(EXT) QNFE 2822 3166 1716 
2 
SUMT CBI .84527E+05 .84526E+05 .84525E+05 
(LIN- NAC Q* Q* O* 
EXT- NFE 578 390 388 
INT) 
3 
SUMIr OBJ . 84480E+05 .84480E+05 .84480E+05 
(QUAD~ NAC O* Ox Ox 
EXT NFE 714 Fie iD 
INT) 
4 
SUMT OBJ .84512E+05 ~84512E+05 .84512E+05 
(CUBIC NAC Q* Q* O* 
EXT NFE 535 535 585 
INT) 
5 OBI .46665E+05 .46849E+05 .43967E+05 
SUMT NAC O* g 2 
(ALM) NFE 4098 3949 5657 


CONMIN RESULTS: 


OBJ: 
NAC: 
NFE: 


popes S405 
~ 
2946 


* = VIOLATED CONSTRAINT(S) 





TABLE XXXVI 


Test Results 234-Bar Space Tower (IONED: 3,8) 


ONE-DIMENSIONAL SEARCH: Golden Section + Cubic Polynomial 
NDV = 56 NCON = 3550 (Stress, Displacement, and Buckling) 
AVERAGE CPU TIME PER FUNCTION EVALUATION: -34968E+00 Seconds 


OPTIMIZER: ie 2 3 4 5 
METHOD OF  FLETCHER- ROBUST D.F.P. Bor.Gas. 
STRATEGY: FEAS. DIR. REEVES M.F.D. 
0 OBJ .58577E+05 .54682E+05 
DIRECT NAC 1 2 
NFE 2608 43 
1 OBJ .51327E+05 .53002E+05 .45573E+05 
SUMT NAC 2 2 O* 
(EXT) NEE 3189 2917 3101 
2 
SUNT OBJ .84529E+05 .84528F+05 .84529E+05 
(LIN- NAC O* O* Q* 
T- NFE 584 585 970 
INT) 
=} 
SUMT OBIT .84480E+05 .84480E+05 .84480E+05 
(QUAD- NAC O* O* O* 
EXT NEE nis 713 713 
INT) 
4 
SUMT OBJ .84512E+05 .84512E+05 .84512E+05 
(CUBIC NAC o* O* o* 
EXT NFE 535 535 535 
INT) 
5 OBJ .46507E+05 .45394E+05 .45364E+05 
SUMT NAC 2* l iL 
(ALM) NFE 4702 4100 4527 


CONMIN RESULTS: 
SBo: .39353E+05 


NAC: 4 
NFE: 2946 


x = JTOLATESMOGONS TRAIN (S) 





TABLE XXXVII 


Test Results 234-Bar Space Tower (IONED: 4,9) 


ONE-DIMENSIONAL SEARCH: 
NDV = 56 NCON = 3550 (Stress, Displacement, and Buckling) 
AVERAG]S CPU TIME PER FUNCTION EVALUATION: ~-32508E+00 Seconds 


Bounds + Polynomial 


OPTIMIZER: i Z 3 4 5 
METHOD OF £FLETCHER- POBUST Der oe BebeGeo. 
STRATEGY : FEAS. DIR. REEVES Mere. Db: 
0 OBJ .51504E+05 -54825E+05 
DIRECT NAC 2 1 
NFE 2279 2070 
a OBJ ~45275E+05 -46030E+05 .35709E+05 
SUM? NAC O* Q* 1* 
(EXT) NEE So 1971 1487 
2 
SUMT OBJ 8452 7E+05 .84526E+05 .84528E+05 
\LIN= NAC O* O* Q* 
=) NEE 529 S56 706 
INT) 
3 
SUMT OBJ -84480E+05 -S4480E+05 .84480E+05 
(QUAD- NAC Ox Ox Ox 
EXT= NEE 698 698 698 
INT) 
4 
SUMT OBJ ~84512E+05 ~84512E+05 .84512E+05 
(CUBIC NAC O* Ox O= 
EXT= NFE 524 De 524 
INT) 
5 OBJ -48 394E+05 -45472E+05 .453218+05 
SUMP NAC 3 L 2 
(ALM) NEE 3423 3769 3598 
CONMIN RESULTS: * = VIOLATED CONSTRAINT (S) 
OBJ : (3935 3E+05 
NAC: 4 
NFE: 2946 
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